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ABSTRACT 



Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular 
' clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The 
' emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are 
detected at 353 GHz and 143 GHz, with amplitudes around -7% and +13%, respectively, indicating that the measured spectra are likely more 
complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 
100 GHz bands, mainly caused by to the contribution of the 7 = 2^1 and 7=1^0 ^^CO and ^^CO emission lines. We derive maps of the 
dust temperature T, the dust spectral emissivity index and the dust optical depth at 250yL/m T250. The temperature map illustrates the cooling 
, of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. 
■ The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant 
' T -j3 anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest 
molecular regions to the faint diffuse regions. We use near-infrared extinction and H i data at 21 -cm to perform a quantitative analysis of the spatial 
variations of the measured dust optical depth at 250 //m per hydrogen atom T25o/A^h- We report an increase of T25o/A^h by a factor of about 2 
between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles. 

Key words. ISM: general, dust, extinction, clouds - Submillimeter: ISM 



1. Introduction 

Planck^ (iTauber et al]l201Ql : lPlanck Collaborat"ionll201 lah is the 
third-generation space mission to measure the anisotropy of the 

^ Planck (http://www.esa.int/P/awc^ ) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries 
Erance and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular resolution from 3' to 5'. The Low Frequency 
Instrument (LEI: iMandolesi et al.ll2QlQl: iBersanelli et al.ll2Ql(i 
iMennella etal.ll2Qllh covers the 30, 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Inst rument (HFI; 
iLamarre et all 120101: iPlanck HFI Core TeamI |201 la) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to O.I K. Polarisation is measured in all but the 
highest two bands (iLeahv et al.1 120101: iRosset et al] l2010l) . A 
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combination of radiative cooling and three mechanical cool- 
ers produces the temperatures needed for the detectors and op- 
tics (Planck Collaboration 2011b). Two data processing centres 
DPCs) check and calibrate the data and make maps of the sky 
Planck HFI Core TeamI [2QTlbl: IZacchei et all l2Qlll) . Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for Galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration, 201 Ih-u. This paper presents the early re- 
sults of the analysis of Planck/HFI observations of the Taurus 
molecular cloud. 

IRAS provided a complete view of the interstellar matter in 
our Galaxy in four photometric bands from 12 to 100//m with 
angular resolution of about 4'. The 12, 25, and 60yum bands 
are mainly sensitive to the emission of the smallest interstel- 
lar grains. They are aromatic hydrocarbons (large molecules) 
and amorphous hydrocarbon grains that undergo stochastic heat- 
ing upon photon absorption, and tend to emit most of their 
energy at wavelengths shortward of lOOyum. The large parti- 
cles have dimensions of the order of lOOnm and make up the 
bulk of the dust mass. They are in equilibrium between ther- 
mal emission and absorption of UV and visible photons from 
incident radiation. Only one IRAS band, at lOOyum, is dom- 
inated by the emission of this dust component. DIRBE and 
FIRAS on board COBE produced all-sky maps at longer wave- 
lengths, with angular resolution lower than IRAS (40' and 7°, 
respectively), which allowed the measurement of dust temper- 
atures and spectral indices. The dust temperature was found 
to be on average -17.5 K (with a sp ectral emissivity index 
= 2) in the diffuse atomic medium (iBoulanger et al.l IT996b 
and t o be lower in m olecular clouds with no embedded bright 
stars (iLagache et al.| [l998). Small patches of molecular clouds 
have b een observed i n more detail from the ground by the 
JCMT (I Johnstone & B allv 1999), by balloon-borne experiments 
PRON AOS (Ristorcelli et"an Il998h Arc heops (iDesert et al.1 
l2008ah and BLAST (Netterfield et al."'2009|), and from space by 
Spitzer (e.g., Flasev et al. 2009) and H erschel (e.g., [An dre et al. 
l2QTQl:lMotte et al.ll2010l:lAbergel et aLllioTolilJuvela et al.ll2011l) . 

With Planck, the whole sky emission of thermal dust is 
mapped from the submillimetre to the millimetre range, with 
angular resolution comparable to IRAS. We have complete and 
unbiased surveys of molecular clouds not only in terms of spa- 
tial coverage, but also in terms of spatial content of the data. 
Unlike ground or ballon observations, the maps contain all an- 
gular scales above the diffraction limit without spatial filtering. 
This is extremely important for Galactic science, because the in- 
terstellar matter contains a large range of intricate spatial scales. 
Moreover, the emission is measured with unprecedented signal- 
to-noise ratio down to the faintest parts surrounding the bright 
and dense regions. Up to now, the rotational transitions of car- 
bon monoxide (CO) were the main tracer of the interstellar mat- 
ter in molecular clouds. With Planck thermal dust becomes a 
new tracer. The spectral coverage and the sensitivity allow for 
each line of sight precise measurements of the temperature, of 
the spectral emissivity index, and of the optical depth, indepen- 
dent of the excitation conditions. 

Molecular clouds present a wide range of physical conditions 
(illumination, density, star-forming activity), so are ideal targets 
for studying the emission properties of dust grains and their evo- 
lution. The goal of this paper is to discuss the results of the early 
analysis of HFI maps of the Taurus complex, whic h is one of the 
neare st giant molecular clouds (J = 140 pc, from iKenyon et al.l 
11994 with low-mass star-forming activity. 



First we present the HFI data and the ancillary observations 
used for this paper (Sect.©. Then we describe the emission spec- 
trum of thermal dust measured by HFI and IRAS, and discuss 
the validity of the measured spectra fitting with one single modi- 
fied blackbody (Sect.©. Maps of the dust temperature and of the 
spectral emissivity index are analysed in Sect.lH Then in Sect.[5] 
we discuss the evolution of the dust optical depth/column den- 
sity conversion factor (r/Nu) from the atomic diffuse regions to 
the molecular regions. 

2. Observations 

2.1. HFI data 

We used the DR2 re lease of the HFI ma ps presented in Fig.[T]and 
begin with Healpix (iGorski et al.ll2005h with A^side = 2048 (pixel 
size 1.7). The data processing and calibration are described in 
iPlanck HFI Core T eam (2011b). An important step is the re- 
moval of the CMB through a needlet internal linear combination 
method. 

The syste matic calibration accuracy of H FI is summarized in 
Table 1 (from lPlanck HFI C ore Team' 201 la*). For the two bands 
at 857 and 545 GHz, the gain calibration is performed using 
FIRAS data. The CMB dipole is used at lower frequencies. The 
systematic errors on the gain calibration are 7% for the 857 and 
545 GHz bands, respectively (estimated using the dispersion in 
different regions of the sky), and about 2% for the other bands. 

We estim ated the statistical noise a s follows (see also 
Appendix B of iPlanck Collaborationll201 ItL for details). Two in- 
dependent maps of the sky have been computed by the DPC from 
the first- and second-half ring of each pointing period. Because 
the coverage is identical for these two maps, the standard de- 
viation of their half difference, (Thr, is equal to the standard 
deviation of the average map. For all bands, we compared the 
computed value of (Thr with the standard deviation (Xref of the 
DR2 map computed within a reference 48' x 48' window cen- 
tered at / = 165?43, b = -21?06, chosen in the lowest part of 
the map (white square in the first image of Fig.[T]). Obviously, 
the values of (Xref give only upper limits on the statistical noise 
because of the contributions of CMB residuals, and of cosmic 
infrared background (CIB) and thermal dust fluctuations that in- 
crease with increasing frequencies. Table 1 shows that for the 
100, 143, and 217 GHz bands, (Thr is almost identical to (Xref. 
Therefore we conclude that 

1. the standard deviation of CMB residuals is lower than the 
statistical noise in all bands; 

2. the standard deviation of the CIB anisotropics (CIB A) ap- 
pears lower than the measured statistical noise in the 100, 
143, and 217 GH z bands (which is compatible with the CIB A 
measurements in IPlanck Collaborationll201 fnl) : 

3. the computed values of (Thr give realistic estimates of the 
statistical noise. 

For this early analysis, we used a constant statistical noise in 
each band, taken equal to the values of cthr in Table 1 . 

2.2. Ancillary data 

We combined the Planck maps with IRAS maps at 100/im 
(3000 GHz), using the IRIS (I mproved Reprocessing of the 
IRAS Survey) data computed by Miville -Deschenes & Lagach3 
(2005). The statistical noise of the lOO^m maps is about 
0.06 MJy sr"^ per pixel (pixel size of 1.8'), while the system- 
atic error in the gain calibration from DIRBE is estimated to 
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Fig. 1. IRAS and HFI maps of the Taurus molecular cloud, in MJy sr ^ . The 48' x 48' reference window is seen on the IRAS map at 
100 Jim (3000 GHz). For all maps the average brightness computed within the reference window is subtracted. o 
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Table 1. Calibration accuracy, statistical noise, standard deviation within the reference window (white square in the 3000 GHz panel 
of Fig.[T]), median brightnesses, and average spectrum within the reference window. 



Frequency (GHz) 


100 


143 


217 


353 


545 


857 


3000 


Calibration accuracy 


2% 


2% 


2% 


2% 


7% 


7% 


13.5% 


Statistical noise from (Tur (MJy sr~^) 


0.014 


0.0086 


0.018 


0.038 


0.049 


0.074 




CTref (MJy sr-^ 


0.013 


0.0078 


0.019 


0.050 


0.13 


0.35 


0.37 


Median brightness (MJy sr~^) 


0.042 


0.094 


0.41 


1.67 


6.37 


17.46 


13.9 


Reference spectrum (MJy sr~^) 


0.010 


0.027 


0.11 


0.43 


1.48 


4.1 


4.8 



be 13.5% (from i Miville-Deschenes & Lagachdl2005h . For the 
Taurus molecular cloud, the lOOyum brightness is in the range 
l-20MJy sr"\ so the statistical noise translates into relative er- 
rors in the range 0.03-6%, well below the systematic error on 
the gain. 

The atomic gas was traced on large scales using the H i data 
at 21 cm taken with the Leiden/ Dwingeloo 25 -m tel escope with 
an angular resolution of 36' by lHartmann & BurtonI (Il997ii) . The 
velocity spacing was 1.03kms"^ and the local standard of rest 
(LSR) velocity range was -450 <Vlsr< 400 km s"^. The data 
were corrected for contamination from stray light radiation to 
the 0.07 K sensitivity level ( Hartmann et al. 1996). 

The large- scale survey in the ^^CO (/ = 1 ^ 0) emission 
line taken by IPame et al.l (l200ll) with the Cf A telescope was 
used to trace the molecular gas. The beam width was 8.8' ±0.2'. 
For the observations of the Taurus molecular cloud, the sampling 
distance was equal to 7.5', the channel width is 0.65 km s"^ and 
the channel rms noise was 0.25 K. The data cubes were trans- 
formed into the velocity-integrated intensity of the line (Wco) 
by integrating the velocity range where the CO emission is 
significantly detected using the moment method proposed by 
[Pame et al.i (2001). The statistical noise level of the Wqo map 
is typically 1 .2 K km s"^ 

In order to analyse the central region of the Taurus molecular 
cloud with the ^^CO / = 1 ^ emission lines, we also used the 
FCR AO survey of 98 deg^ conducted with an angu lar resolution 
47" (iNaravanan et al.ll2008l: [Goldsmith et al.ll2008l). We appHed 
the error beam scaling factor recently proposed by Pineda et al. 
(I2Q1Q) . so that the intensities of the Cf A and FCRAO surveys of 
the ^^CO emission line are compatible. The statistical noise per 
pixel (with pixel size of 0.33') of the velocity-integrated inten- 
sity maps is about 0.4 K km s"^ for the ^^CO line. 

The column density can also be traced from the near-infrared 
(NIR) extinction using the 2MASS point source catalogue (e.g., 
iDobashi et al.l l2QQ5l: iPineda et al.l [2010). For comparison with 
HFI data, we used the extinction map of the ce ntral part of 
the Ta urus molecular cloud recently created by [Pineda et"aD 
([20TqI) and shown in Fig.O It is Nyquist- sampled with an an- 
gular resolution of 200", and corrected for the contribution of 
at omic gas to the total extincti on. Using the extinction curve 
of [Weingartner & Draind ([200 1[) and adopting the standard ra- 
tio of selective to total extinction Rv=3.l for the diff'use ISM 
([Savage & Mathis 1979), the infrared colours are converted to 
visible extinction Ay or column densities A^h (= Ay x 1 .87 x 10^^ 
cm^). Because the number of background stars used to compute 
the extinction decreases with increasing extinction, the error in 
Ay increases from about 0.2 mag at low extinctions (Ay = 0- 
Imag) to 0.5 mag at high extinctions (Ay ~ 10 mag), with an 
average error of 0.29 mag. For higher extinctions (Ay > 10 mag) 
the extinction map only gives lower limits. 
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Fig. 3. Spectrum of one pixel in the non-molecular region: 
crosses, total brightness /tot; triangles, average brightness /ref 
within the reference 48' x 48' window (seen on the 3000 GHz 
panel of Fig. [T]); squares, /tot - /ref; and diamonds, CIB spectrum 
at th e central frequency of the HFI filters from 857 to 217 GHz, 
from lFixsen etaLl ([1998[) . 



3. Emission spectrum of thermal dust 

3.1. Reference spectrum 

The combination of IRAS and Planck data provides the spectral 
energy distribution (SED) from 3000 GHz (lOOyum) to 100 GHz 
(3 mm) for each pixel of the maps. The CMB has been removed 
from the maps we use. The maps could contain some CMB resid- 
ual, but with an amplitude lower than the statistical noise even 
in the low frequency channels (see Sect. l2.1b . 

In all bands, the maps contain Galactic and non-Galactic 
emission that is not associated with the Taurus complex. 
Therefore, we subtracted for all maps the average brightness 
computed within the reference 48' x 48' window, chosen in the 
faintest region (white square in the first image of Fig. [T] and cen- 
tral coordinates given above in Sect. l2.1b of the map. This is il- 
lustrated for one pixel in Fig.O The average spectrum in the 
reference window is also shown in Fig.[3l together with the CIB 
spectrum (from FIRAS data by Fixsen et al. 1998) at the central 
frequencies of the HFI filters. We see that the reference spec- 
trum is mainly caused by Galactic atomic emission (there is no 
detected emission within the ^^CO / = 1 ^ line), since the 
CIB contributes not more than 5-10%. 

In order to derive the dust optical depth per column density 
in the atomic phase (Sect. l5.2l) , the emission measured within the 
same reference window will be subtracted from the H i data. 
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Fig. 2. Upper left panel: column density derived from the Hi data at 21 cm (iHartmann & Burtonll 19971) . Upper right panel: ^^CO 
(/ = 1^0) velocity integrated emission line (Dame et al. 2001). Lower left panel: ^^CO / = 1 ^ velocity integrated emission 
line (iNaravanan et al.ll2008l:lGoldsmith et al.ll2008 ). Lower right panel: NIR extinction using 2MASS (iPineda et al .11201 Ol) . The ^^CO 
and NIR extinction maps are smoothed at the angular resolution of the 143 GHz HFI band (FWHM: 7.08'). 



3.2. Choice of the spectral bands 

We focus our analysis on the emission of dust particles in equi- 
librium between thermal emission and absorption of UV and vis- 
ible photons from incident radiation. Therefore, we did not use 
the IRAS maps at 12, 25, and 60yum because of the contribution 
of small dust particles transiently heated each time they absorb 
a UV/visible photon. In all spectra presented in this paper, we 
left the 60 jim data points in the figures to illustrate this contri- 
bution, but these data points were not used to analyse the SEDs. 
In the lOOyum band, the contribution of these small particles is 
expected to be lower than 10% if the intensity of the interstellar 
radiation field is of the order of the value in the local diff'use ISM 
(e.g., Compiegneet al. 201 1), as is the case in the Taurus molec- 
ular complex. This contribution is in any case lower than the gain 
calibration error of IRAS at 100 yum (13.5%, see Sect.O. 

For this early analysis, we did not use the data taken at 100 
and 217 GHz to analyse the thermal dust emission because they 
are expected to be contaminated in the molecular regions by ro- 
tational / = 1 ^ and / = 2 ^ 1 ^^CO emission lines and 
by ^^CO emission lines with a lower amplitude. Higher /-lines 
in the other bands have a lower relative amplitude and are ne- 
glected for this early analysis. We also neglected the emission of 
molecular lines tracing the densest regions. Therefore, in a first 



step only the bands at 3000, 857, 545, 353, and 143 GHz (100, 
350, 545, 850, and 2100 yum) were used to analyse the thermal 
dust emission. 



In practice, we smoothed all maps to the angular res- 
olution of the 143 GHz data (FWHM: 7.08'), assuming 
isotropic Gaussian beams a n d tak ing the FWHM from 
Miville-Deschenes & Lagache (120051) for IRAS and from 
Planc k HFI Core Team (20r\J) for HFI. We did not take into ac - 
count the ellipticity of the PSF (iPlanck HFI Core Teamll2011ah . 
This may introduce some error for the detailed analysis of indi- 
vidual structures, but our goal is not to to derive any quantita- 
tive results for individual structures, but to extract from a pixel 
by pixel analysis some quantitative information on the emission 
that emerges from the diff'erent phases. The 100 GHz map will 
be used to compute the residuals to the fits (see Sect. l3.5l) , and is 
left at its original resolution (9.37). Figure |4] shows the spectra 
taken at positions in the non-molecular region and in the molec- 
ular region. 
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Fig. 4. SED of two pixels (top: non-molecular region, bottom: 
molecular region). The squares are data, the solid line is the fitted 
model, and the crosses are the fitted model integrated within the 
bands. The fits are performed using the 100 yum, 857, 545, 353 
and 143 GHz bands (red squares), and using the statistical noise 
discussed in Sect.O which is too low to be visible on the figure. 
Significant excess in the 217 and 100 GHz bands caused by ^^CO 
and ^^CO emissions are detected in the molecular spectrum. The 
60 Jim data points are not used to analyse the SED because of the 
contribution of small dust particles transiently heated each time 
they absorb a UV/visible photon. 

3.3. Principle of fitting 

In this early analysis, the fitting function is a single and optically 
thin modified blackbody: 

Iv = T,,x|^j^x W), (1) 

where Ty^ is the dust optical depth at frequency vo, is the spec- 
tral emissivity index. By is the Planck function, and T the dust 
temperature. All quantitative values of the dust optical depth will 
be given at the frequency vq = 1200 GHz (250yum) to be compa- 
rable to previous analyses, and we define Ty^ = T250. 

The three computed parameters for each pixel are T, JS, and 
T250- We used the IDL MPFIT routine, w hich performs wei ghted 
least-squares curve fitting of the data (iMarkwardtl l2009b tak- 
ing into account the noise (statistical noise or calibration un- 
certainty) for each spectral band. We applied colour-correction 
factors computed using version 1 of the transmission curves 
(iPlanck HFI Core Teaml[20TTd) . 



The goal of this adjustment is to reduce the properties of the 
SED measured for each line of sight to the set of three param- 
eters T,j3, and T250. The temperature, spectral emissivity index, 
and optical depth of the emitting dust particles along the line of 
sight can obviously vary, so the fitted values of the three param- 
eters, while representing some average properties, cannot give 
the complete picture of the dust particles located along that line 
of sight. Moreover, we assumed that the spectral emissivity in- 
dex of the measured spectra is constant from far-infrared (FIR) 
to millimetre wavelengths. The spectral emissivity index of dust 
emission may vary due to temperature-dependent mechanisms 
at low temperatures, including free-charge carrier processes, 
two-phonon diff^erenc e processes, and absorption mechanisni s 
in two-level systems (lAgladze et"aDll996l: iMenv et al.1 l2007h . 
Moreover, the measured spectra can be broadened around the 
peak of the modified blackbody in the submillimeter (submm) 
because of the contribution of dust at diff'erent temperatures. 
Thus, the spectral emissivity index j3 of the measured spectra 
co uld increase from th e submm to the millimeter, as illustrated 
bv lShettv et aP (l2009bh . 

3.4. Example of spectra. Systematic errors on the 
parameters derived from the fit 

Figure|4]is an illustration of the fitting for two pixels, one taken at 
a position with detected CO emission (within the ^^CO / = 1 ^ 
line), the other at a position with no detected CO emission. We 
used for the fitting the statistical noise for HFI and IRAS data 
discussed in Sect.O taking into account the smoothing of the 
map at the angular resolution of the 143 GHz band. 

Simulations have been performed to understand and to quan- 
tify the propagation of the calibration errors, the statistical 
noise, and the CIBA in the determination of and T250 (see 
AppendixlAl for details). The calibration errors propagate into 
systematic errors in T,/3, and T250 of about 0.7 K, 0.07, and 18%, 
respectively. We have also shown that the statistical noise and 
the CIBA propagates into statistical noise levels in T,j3, and T250 
in the range 0.1-1 K, 0.025-0.25 and 2-20%, respectively, de- 
pending on the 100 yum brightness (10-1 MJy sr"^). Moreover, 
both for systematic errors and statistical noises, the three param- 
eters are strongly correlated or anti-correlated (Fig. lA.lb . 

The two spectra of Fig.|4l show that to first order, a sin- 
gle modified blackbody gives a reasonable representation of the 
SED measured by IRAS and HFI. Obviously by increasing the 
number of free parameters in our fit (e.g., by using two modified 
blackbodies with diff'erent temperatures and spectral indices) it 
is possible to improve the fits significantly, but this is not our 
goal in this early paper. We used exactly the same method to fit 
the spectra for all pixels of the map to derive the temperature 
map and the spectral emissivity index map shown in Fig.[7l 

3.5. Analysis of the residuals 

The fit residuals allow us to assess the limitations of usi ng a 
single modified blackbody to fit the data. As seen in Sect. l3.2l 
fits were made to the 3000, 857, 545, 353, and 143 GHz data 
all smoothed to the angular resolution of 143 GHz. The residual 
map at 100 GHz was computed from the diff'erence between the 
synthetic spectra smoothed at the 100 GHz resolution and the 
data at 100 GHz. The residual map of the fitting in all bands is 
shown in Fig. [51 Figure[6] shows the relative residual map (resid- 
ual map divided by measured map). 
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3.5.1 . Residuals for the five fitted bands 

We see in Fig. [6] that the relative residuals are distributed around 
zero and below 1-3% for the 3000, 857, and 545 Hz bands. On 
the other hand, the residuals at 353 GHz are systematically neg- 
ative, with a median relative amplitude about -7%, while the 
residuals at 143 GHz are systematically positive, with a median 
relative amplitude about +13%. As can be seen in Figs.[5]and[6l 
these residuals are spatially correlated with the measured bright- 
ness, so they are related to the emission spectrum of the dust 
particles located in the complex. 

We have seen in Sect. l2.1l (Table 1, from 
iPlanckHFI CoreTeamI 1201 lah that the calibration errors 
are estimated to be 13.5% at 3000 GHz (100 yum), 7% at 857 and 
545 GHz, and 2% at 353 and 143 GHz. These numbers translate 
into systematic errors of about 5% and 4% on the 353 and 
143 GHz brightnesses (computed from the fit of the five bands), 
respectively. We conclude in this early analysis that the negative 
residuals at 353 GHz (with a median amplitude about -7%) 
and the positive residuals at 143 GHz (with a median amplitude 
about +13%) are not caused by calibration errors. The negative 
residuals at 353 GHz could be caused by a broadening of the 
measured spectra around the peak of the modified blackbody 
(3000-545 GHz spectral range) because of the contribution of 
dust at different temperatures along the line of sight. On the 
other hand, the positive residuals at 143 GHz suggest a flattening 
of the emission spectra at low frequencies. 

3.5.2. Residuals at 21 7 and 1 00 GHz 

As shown in Fig. [51 there is a striking spatial correlation between 
the residual maps at 100 GHz and 217 GHz and the map of the 
integrated emission of the ^^CO / = 1 ^ emission line. This 
confirms that 100 GHz and 217 GHz residuals are dominated by 
a contamination by CO molecular lines (mainly ^^CO, but also 
^^CO and other isotopes and molecules in the densest regions). 
However, we have seen above that the measured dust spectra are 
more complex than a single modified blackbody from 3000 GHz 
to 100 GHz, so the residuals computed at 100 GHz and 217 GHz 
in this early analysis and shown in Fig. [5] are only indicative of 
the CO emission. 

4. Analysis of the dust temperature and spectral 
emissivity index maps. 

4.1. Dust temperature map 

The dust temperature map is shown in Fig.O The black regions 
are artifacts: they correspond to low values of the temperature 
(below 12K) caused by statistical noise and CIBA in the faint 
regions (with /(lOOyum) < 1 MJy sr"^ as also illustrated in our 
simulations presented in Fig. lA.2l) . Excluding these regions and 
some residual stripping in the faintest regions with amplitudes 
around 0.15 K, most of the variations seen on the temperature 
map are real because their amplitude is higher than the noise 
computed from our simulations (about 0.1-1 K for / (lOOyum) = 
10-1 MJy sr-^ from ¥ig\K^. 

At least three regions with diff'erent temperature distribu- 
tions can be identified in the temperature map (Fig.jTl) and its 
histogram (Fig. [8]): 

- the outer parts of the molecular cloud (with no detected ^^CO 
emission), with temperatures ~ 16-17.5 K; 

- the molecular cloud with detected ^^CO emission but no de- 
tected ^^CO emission, with temperatures ~ 15-17 K; 




^A\ ^ ^ I 
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Fig. 9. Correlation between the spectral emissivity index and 
the dust temperature. Upper panel: for all pixels. Lower panel: 
for pixels in the molecular region with detected CO emission 
and/(100yum) > lOMJysr"^ The sol id and dashed lines show 
the relations deduced from Archeops (D esert et"aD l2008bl) and 
PRONAOS (Pupae et al 20031 respectively. 



- the densest parts of the molecular cloud with detected ^^CO 
emission which coincide with the well-known dense fila- 
ments, with temperatures ~ 13-16 K. 



Most of the diff'erences between the IRAS map at 100 yum and the 
HFI maps (Fig.[T]) are caused by temperature variations. Except 
in diff'use regions with uniform illumination and constant dust 
temperature, this demonstrates that IRAS maps at lOOyum used 
alone cannot properly probe the gas column density. Comparable 
temperature variations have also been found by iFlagey et"aD 
( 2009) in the central region of our field (average temperature 
about 14.5 K with a dispersion of 1 K) by combining Spitzer 
160yum and IRAS lOOyum maps and assuming a constant value 
of p equal to 2. 
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Fig. 7. Left panel: dust temperature map. Right panel: spectral emissivity index map. 





Fig. 8. Dust temperature and spectral emissivity index histograms: black, all pixels; green, pixels without detected ^^CO emis sion; 
red, pixels with detected ^^CO emission but no detected ^^CO emission in the central molecular region covered by .Pineda et al.l 
(I2Q1QI) : and blue, pixels in the central molecular region with detected ^^CO emission. 



4.2. Spectral emissivity index map 

The spectral emissivity index presents a symmetric distribution 
(Fig. [8]) with average and median values both equal to 1.78 and a 
standard deviation cr^ = 0.08. The standard deviation caused by 
statistical noise and CIBA estimated from our simulations is in 
the range 0.025-0.25 for 7 (100 yum) = 10- 1 MJy sr'^ (Fig|A2]). 
Moreover, we can see in Fig. [7] and more explicitly in the first 
panel of Fig. [9] that there is a some anti-correlation between T 
and 13. The higher values of (> 1.8) are found in the coldest 
(about 14 K) structures. 

We sho w in A ppendixlAl and it was also discussed in de- 
tail by S hetty et al.l (l2009 a,b) that the instrumental errors always 
produce intrinsic anti-correlation between T and JS. A significant 
fraction of the anti-correlation seen in Fig. [9] could be caused by 
noise. However, high values of j3 generally correspond to bright 
regions. This is illustrated in the second panel of Fig.[9l which 
shows the T - JS correlation diagram for pixels in the molecular 
region with detected CO emission and 7 (100 jim) > 10 MJy sr"^ : 
the r - anti-correlation is still visible, with an amplitude sig- 
nificantly higher to what is expected from the contribution of the 
statistical noise and the CIBA computed from simulations (green 
points in the first panel of Fig. lA.lK 



Previous observations of thermal dust emission at FIR to 
millimetre wavelengths in a variety of Galactic regions indi- 
cate an anti-correlation between the fitted values of the spectral 
emissivity index jS and the dust temper ature, f rom PRONAOS 
data for T = 12 - 20 K (Dupac e tjpl2003l) . Archeops data 
for r = 7 - 27 K (Desert et al. 2008bX and Herschel data fo r 
r = 10 - 30K (e.g., Anderson et al. 2010; Pa radis et al.ll2010h . 
We see in Fig.[9]that the T -yS anti-correlation we find is steeper 
than the relation deduced from PRONAOS, but comparable to 
the relation deduced from Archeops. This T-p anti-correlation is 
for line-of-sight-averaged data and reveals an intrin sic property 
of the emission spectrum of the dust particles (e.g. jMeny et al.l 
l2Q07l:lBoudet et al.ll2005l) . The amplitude of the anti-correlation 
may be higher for the emission spectrum of the dust particles 
than for line-of-sight-averaged data because in the general case 
there is a combination of dust temperatures in each line of sight. 



4.3. Comparison with other Planck results and dust models 



As shown in iPlanck Collabora tion (20111), both FIRAS and 
UFl+IRAS spectra of the local diff'use ISM can be well-fit by 
a modified blackbody with T = 17.9 K and j3 = 1.8. Morever, 
median values of T = 17.7 K and JS = 1.8 are also found at 
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b > 10 ° in the all-sky analysis usin^ the 3000, 857, and 545 GHz 
bands (iPlanck Col laboration 201 lo). These two results are fully 
compatible with our distribution of spectral emissivity indices 
(Fig. [8]) centred on 1.78 with a standard deviation of 0.08. We 
also saw in Sect. l3.4l that the systematic error in (3 is estimated to 
be 0.07. 

We compa re here oiir results to the dust mod- 
els of .Draine & Lil (l2007l) and ICompiegne et al.l 
(l2QTlh . which were c omputed with the DustEM tool 
(|http://www.ias.u-psud.fr/DUSTEM). As far as the proper- 
ties of thermal dust are concerned, silicates are the same 
in both models, whereas the carbon grains are taken to be 
grap hite (Draine & Li 2007) , or hydrogenated amorphous car- 
bon (ICompiegne et al.ll201 ll) . We first computed dust emission 
models for all grain types in the diff'use ISM, i.e., hea ted by 
the sta ndard interstellar radiation field (ISRF) of Mathis etaH 
([l983h . Next, we generated simulated data points by taking the 
flux densities of our model spectra in the bands at 100 yum, 857, 
545, 353 and 143 GHz. We then applied our fitting procedure 
to these simiilated band flux densities. Fitting the model of 
iDiaine&Li (l2007h yields T ^ 19.6 K and ^ 1.67. The 
latter value is intermediate between that obtained from fits for 
silicates (1.57) and grap hite (1.93) alone. Fitting the model of 
ICompiegne et al.1 (120111) yields T ^ 20.5 K and /? ^ 1.52 (see 
Fig.[TO|). Individual fits for silicates and amorphous carbon give 
similar values of p. We note that the highest residuals from a 
single modified blackbody are found between 150 and 300 yum, 

in the spectral gap between IRAS and H FI (Fi g [TOb. 

For both the iDraine & lH (l2007h and ICompiegne et al.1 

(I2OIII) models, the fitted values of p are slightly below the cen- 
tral value found in the Taurus molecular cloud (1.78, with a sys- 
tematical error of 0.07). However, it is worth noting that to ac- 
count for the FIRAS spectrum of the diff'use ISM, Li & Draind 
(I2OOII) decreased p from 2.1 to 1.6 in the absorption efficiency 
of silicates at wavelengths above 250yum. If we do not ap- 
ply this correction and simply extrapolate the absorption ef- 
ficiency of silicates as a single power law with p - 2.1 for 
/I > 30 yum, our fit provides T ^ 1 9.4 K and p ^ 1.73 for 



the model of ICompiegne et al.1 (1201 ll). and T ^ 18.2K and 
P ^ 1.97 for the model of iDraine & Li (2007). The value p 
obtained with the ICompiegne et al.. ( 2011) model is then com- 
patible with th e observations. Finally, th e dust temperature de- 
rived from the 'Compieg ne et al.l (1201 ll) model is sig nificantly 
highe r than the value of 17.9 K found by iPlanck Co llaboration 
(l2Qlltb in diffuse cl ouds: this is because the amorphous carbon 
of ICom piegne et al.] f2011) absorbs the ISRF in the near-IR (be- 
tween 1 and lOyum) more efliiciently than graphite or silicates. 

This first comparison of the Planck data and recent dust mod- 
els shows that the emission of thermal dust can be represented 
as a first approximation by a single modified blackbody with a 
constant spectral emissivity index p. 



5. Dust optical depth per unit column density 

We have seen in the previous section that the spectral cover- 
age of Planck/HFl allows for the measurement of the dust tem- 
perature and of the spectral emissivity index p for each line of 
sight. The third adjustable parameter is the dust optical depth at 
250 yum, T250. To allow a direct comparison of dust optical depth 
results with the other Planck Early Results papers, we fitted the 
data using only the bands at 3000, 857, and 545 GHz, holding p 
fixed at 1.8. Compared to fits with five bands (including 353 and 
143 GHz), the fitted optical depths do not change by more than 
about 10%. 



For the first time we have an unbiased map of the optical 
depth of thermal dust in a molecular complex from the most dif- 
fuse regions to the densest parts. The optical depth map in Fig.fTTI 
has a dynamic range greater than 100. The statistical error and 
the systematic noise in T250 are estimated to be about 1-10% and 
12%, respectively, using the method described in AppendixlAl 
(but with a fixed value ofp= 1.8, and considering fitting of only 
the three bands at 3000, 857, 545 GHz). 

5.1. Independent tracers of the column density 
5.1.1. Atomic phase 

The column density of the atomic gas is traced using the Hi 
data at 21 cm taken with the Leide n/Dwingeloo 25 -m telescop e 
with an angular resolution of 36' (IHartmann & Burtonlll997l) . 
In the optically thin hypothesis, the velocity-integrated emis- 
sion can be converted to column density using the classical fac- 
tor lK"^km"^s = 1.81 X 10^^ cm"^. However, the Hi emis- 
sion is subject to self- absorption in the cold neutral medium 
(CNM), and Nu can be underestimated. An exact calculation 
requires knowledge of the density profiles of the CNM com- 
ponents alo ng each line of si ght. This estimate has been per- 
formed by iHeiles & Trolandl (120031) by measuring the emis- 
sion/absorption of the 21 -cm line against a number of con- 
tinuum sources. Correction factors around 1.25 are found in 
the Taurus/Perseus region. In our case, the precise correction 
pixel by pixel for each line of sight is n ot possible. We tested 
the simple correction method discussed in lPlanck CoUaborationl 
(201 It), which assumes a constant spin temperature Ts (this is 
not really justified because Ts varies from the warm neutral me- 
dia (WNM) to the CNM). Correction factors of about 1.1-1.5 
are obtained for Ts = 100 K. In this early analysis we decided to 
apply a constant multiplicative correction factor of 1.25, with a 
conservative uncertainty of 20%, to derive from the H i data the 
column density map for the atomic phase (Fig. [21). In any case, 
we checked that a diff'erent choice for the correction method does 
not aff'ect the quantitative analysis of the dust optical depth per 
unit column density in the molecular phase (see below). 

The Taurus molecular cloud is 15° from the Galactic Plane 
and contains large-scale emission associated with background 
atomic gas with LSR velocities from -50kms"^ to Okms"^. 
This explains the north- south Galactic gradient in the HFI and 
IRAS maps (Fig.[T]) and also in the dust optical depth map 
(Fig.[TT]). An east- west gradient caused by a filamentary struc- 
ture that extends away from the Galactic plane and crosses the 
eastern part of Taurus is also detected, with velocity in the same 
range as the velocity in the ^^C O / = 1 ^ e mission line, i.e., 
from ~ Okms"^ to ~ 15kms"^ (iNarayanan et al 2008.) . 



5.1.2. Molecular phase 

The column density associated with the molecula r regions can 
be traced with the NIR extinction map of Pineda et al.l |201Q| 
(Fig. 121), which gives the extinction caused by dust associ- 
ated with th e non-atomic component. Special care is taken by 
Pineda et al. 1 (f2010) to remove the overall extinction associated 
with H I located between the background stars used and the Earth 
(0.3 mag), and also the extinction associated with the widespread 
Hi emission (0.12 mag). 

The / = 1 ^ line of ^^CO is a widely-used tracer of 
the molecular phase. However, this line is sensitive to variations 
in abundance (depletion in densest regions, formation and de- 
struction), excitation conditions, and radiative transfer eff'ects 
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Fig. 10. Left panel: fit of the DustEM mod el of ICompiegne et alJ (l201lh for the diff'use ISM heated by the standard interstellar 
radiation field (ISRF) of iMathis et all (Il983l) . The black solid line is the model, the squares are the model in the photometric bands 
at 100 yum, 857, 545 and 143 GHz, and the red solid line is the fitted spectrum. Right panel: relative residuals of the fit. The solid 
line is the conti nuous model and th e triangles are the model in the photometric bands. The dashed line shows the relative residuals 
of the fit for the Draine & 11 (120071) model. The increase of the residuals at wavelengths below lOOyum is caused by the contribution 
of transiently heated small particles. 
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Fig. 11. Maps of the dust optical depth at 250yum (1200 GHz). Left panel: total optical depth derived from the pixel-by-pixel fit of 
the HFI and IRAS data. Right pan el : o ptical depth map of the molecular phase alone (the optical depth associated with the atomic 
phase has been removed, see Sect. l5.3l) . 



(the line is generally optically thick), which explains the dif- 
feren ce between maps of integrated ^^C0(/ =1^0) emis- 
sion (i Dame et al.ll200H Fig.[2]) and our dust optical depth map 
(Fig.[TT]). A detailed comparison between the NIR extinction and 
the integrate d emission of the ^^ CO (/ =1^0) line has been 
presented by lPineda et al.l (|201Q|). In this paper, our strategy is to 
use the NIR extinction map of iPineda et al. l (l20T0h as a quanti- 
tative tracer of the column density of the molecular phase. We 
also used th e ^^CO (/ = 1 ^ 0) integrated emission map of 
iDame et aL 1 (^01), which covers the whole complex, to define 
regions containing (or not containing) molecular material. 

5.2. Dust optical depth per unit column density in the atomic 
medium 

To measure the dust optical depth at 250yum per unit column 
density T250 /Nu in the atomic phase, we used the column den- 



sity map derived from H i data presented in Fig.[2l We subtracted 
from this map the averaged column density computed within 
the reference window (see Sect. l3.1b . Then we smoothed our 
dust optical depth map (Fig.[TT]) to the angular resolution of the 
Hi data (FWHM 36') assuming Gaussian beams. We also se- 
lected all pixels with no detected CO emission, using the ^^CO 
/ = 1 ^ map smoothed to the H i resolution and with the 
criterion Wqo < 0-5 Kkms"^ 

The dust optical depth at 250yum T250 as a function of the 
atomic column density for the selected pixels is shown in Fig.[T2l 
The relationship presents some dispersion, which may be owing 
to the uncertainti es in t he Hi opacity correction (estimated to 
be 20 %, see Sect. l5.1.ll , to the statistical noise in the computed 
values of T250 (about 1-10%, see above), and also to the contri- 
bution of some molecular material not detected on the ^^CO sur- 
vey. However, the relationship appears linear over the full range 
of column density from 1 x 10^^ cm"^ to 3 x 10^^ cm"^, with a 
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Fig. 12. Dust optical depth at 250 jim as a function of the atomic 
column density, for pixels with no detected CO emission. The 
red line shows the result of the linear regression: T250 = 1.14 x 
10-25 xA^H + 4.1 xlO-^ 
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Fig. 13. Dust optical depth at 250yum as a function of the col- 
umn density A^h computed from the NIR extinction map of 
iPineda et"aDl2010 ' (shown on the lower right panel of our Fig. [2]), 
for pixels with detected CO emission {W (CO) > SKkms"^). 
The black line shows the result of the linear regression: T250 = 
2.32 X 10-25 X A^H - 1.44 X lO'^. 



5.3. Dust optical depth per unit column density in the 
nnolecular phase 

In order to measure T250 /A^h in the molecular phase, we need to 
subtract the dust optical depth associated with the atomic phase 
from the dust optical depth map. In this paper the atomic gas is 
traced using Hi data, which have an angular resolution of 36', 
which is significantly lower than Planck/HFl. However, the lo- 
cal densities are higher in the molecular phase than in the atomic 
one, so we can consider that in lines of sight containing detected 
CO emission the small scale spatial fluctuations in the dust op- 
tical depth maps are dominated by density fluctuations in the 
molecular phase. Consequently, first we computed the optical 
depth map of the dust associated to the atomic medium alone 
T25o,Hi by multiplying the Hi column density by the value of 
found in the previous section. Then we subtracted the 
T25o,Hi map from the dust optical depth map to obtain the map 
of the optical depth of the dust associated to the molecular phase 
(Fig. [11]). We see that the north- south and east- west gradients de- 
tectable in the total dust optical depth map have disappeared, as 
expected. 

The dust optical depth at 250 fim T250 in the molecular phase 
as a function of th e column density A ^h computed from the NIR 
extinction map of I Pineda et al.ll201Q (shown in the lower right 
panel of our Fig. [2]) is presented in Fig. [131 The statistical and 
systematic errors on T250 are estimated to be about 1-10% and 
12%, respectively (see Sect.O. The statistical error on A^h is in 
the range 0.2 - 0.5 x 10^^ cm^ for A^h = - 20 x 10^^ cm^ (see 
Sect. 12.21 and taking into account the smoothing at the angu- 
lar resolution of the 143 GHz HFI band). The systematic error 
on A^H is additive and is mainly caused by the uncertainty on 
the zero level from the correction of the extinction associated to 
the atomic phase located between the background stars and the 
Earth , and is estimated to be 0.5 - 1 x 10^^ cm^ by lPineda et al.l 
I2OIOI 

We see in Fig. [13] that the T250 - A^h relationship appears lin- 
ear. The dispersion is mainly caused by the statistical errors on 
A^H (0.2 - 0.5 X 10^^ cm^) and T250 (1-10%). The slope gives a 
measurement of the dust optical depth per unit column density 
in the molecular phase T25o/A^h = 2.32 ± 0.3 x 10"^^. The uncer- 
tainty takes into account the statistical error on T250 and A^h, and 
the systematic error on T250 (12%). 

The T25o-A^H Hnear relationship presents a faint non-zero pos- 
itive residual A^h (^250 = 0) ^ 6 x 10^^ cm^ compatible with the 
uncertainty on the zero level of the A^h niap. 

5.4. Discussion 

5.4.1 . Previous observations 



slope of T25o/A^H = 1.14 ± 0.2 X 10"^^ cm^. The uncertainty of 
'^250 /Nu takes into account the uncertainty in the opacity correc- 
tion of H I data, the statistical noise and the systematic error in 
the dust optical depth map (12%). 

We conclude that the value of T250/A/H computed in the 
atomic medium in the Taurus region appears to be consistent 
with the standard valu e for the diff'use ISM, 1 x 10"^^ cm^ 
(iBoulanger et al.lll996h . 



The possibility of a higher value of T25o/A^h in dense clouds 
than in the diff'use ISM has been a long-stan ding question 
(lOssenkopf & Henningll994lHenning et al This was not 

unexpected, as dust grains must go through coagulation pro- 
cesses in dense environments, which tends to increase r/Nu 
(e.g., Ossenkopf & Henning 1994; Stognienko et al. 1995). Low 
grain temperatures (~ 13K) have been observed also with the 
PRONAOS b alloon in relatively d ifi'use clouds, the translucent 
Polaris Flare ([Bernard et al.lll999b and on e molecular filament 
(for Ay > 2) in the Taurus molecular cloud (Stepnik et al. 2003). 
As demonstrated by the authors, this low temperature cannot be 
explained by the eff'ect of extinction alone: one needs to increase 
the value of T250 /A^h (with A^h derived from the NIR extinction) 
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by a factor of about 3 compared to the standard value for the 
diffuse ISM of 1 x lO'^^ cm^. 

An increase of T25o/A^h in the FIR by a factor > 1.5- 
4 (always with A^h derived from the NIR extinction) has 
been detected from ISQP HOT ob servations of s everal high- 
latitu d e translucent clo i ids ([Cambresv et al. 2 Q0T1: iBurgo et al 



20031: iRidderstad et alJ l2006l: iKiss et alJ l2Q06iriLehtinen et al , 
20071) and TMC-2 d Burgo & Laureiisl l2005'). b y Spitzer o bser- 



vations of the Perseus molecular cloud (Schnee et al .1120081) and 
of the Taurus molecular cloud (Flagey et alj 120091) . and more 
rencentely by Her schel observations of Galactic dense cores 
(lJuvela etal.lf201lh . Because the decrease in temperature is gen- 
erally observed to be associated with a decrease of the 60 yum 
over lOOyum intensity ratio, which traces the abundance of small 
grains relative to that of big grains, coagulation has often been 
invoked to e xplain this emissivity in crease (Bernard et al.ll 19991: 
IStepnik et al.]l2003l: ICambresv et al.ll2005h . 

However, the increase of T25o/Nu with increasing column 
density is not systematically observed. No excess is detected 
in the Corona Australis molecular cloud from Spitzer/MIFS or 
APEX/Laboca data using als o the NIR e xtinction as a tracer 
of the column density ( Juvela "etaD l2009l) . Finally, contradic- 
tory results have also been obtained using the ^^CO(J = 
1 0) emission line to t race t he column density. For in- 
stance, llRoman-Duval et all (l2010l) found that the deviations 
between CO surface density and FIR emission (measured by 
Herschel/SFIRE and Spitzer/MIFS) are more probably caused 
by H2 envelopes not traced by CO and therefore not accounted 
in the column density, th an by gas-to-dust ra tio or T25o/A^h 
variations. On large scales, iParadis et al.l (l2009l) conclude from 
DIRBE, Archeops, and WMAP data that the dust that gives an 
excess of T25o/A^h in the submm may recover its diffuse ISM 
value in the millimetre. This contradicts current scenarios of dust 
coagulation, for which a constant emissivity increase over the 
whole FIR- submm wavelength range is predicted for aggregates 
of astronomical silicate grains. Dust optical constants in the FIR- 
submillimetre range, however, are not well constrained by labo- 
ratory studies. 

In any case, instruments before Planck never had the angular 
resolution, appropriate spectral coverage, sensitivity, and map- 
ping capability to perform full and unbiased surveys of the ther- 
mal dust emission within individual complexes, from the most 
diffuse regions to the densest parts. The key observational ques- 
tions are: where and on which angular scale do the dust proper- 
ties evolve in the ISM. Resolving these questions is the first step 
in understanding the physical processes that regulate the optical 
properties of thermal dust. 



5.4.2. Planck results 

We have seen in Sections [S!2l and [53l that the averaged value of 
the dust optical depth at 250 yum per unit column density T250 /Nu 
increases from 1.14 ± 0.2 x 10"^^ cm^ in the atomic phase to 
2.32 ± 0.3 X 10"^^ cm^ in the molecular phase. It is interesting to 
note that the same value of T250 /A^h was found by iFlagey et aP 
( 2009) in the molecular central region of the Taurus complex 
from IRAS 100 jim and Spitzer 160 jim maps. 

In the molecular phase, the T25o-A^h relationship appears lin- 
ear for A^H = - 15 X 10^1 cm^ (Fig.IH]), so the value of T25o/A^h 
does not appear to depend on the column density A^h- However, 
the value of T250 /Nu for dust particles located in dense regions 
could be higher, due to different effects: 



1. Because there is no embedded heating star, the distribution 
of dust temperatures along a "cold" line of sight, with a 
low value of the measured temperature, should be gener- 
ally broader than along lines of sight with higher measured 
temperatures. Moreover, the measured temperatures for cold 
lines of sight are always warmer than the average tempera- 
ture along that line of sight, so the values of T250 for cold 
lines of sight may un derestimate the average values for dust 
particles B (see also ICambresy et al.l l200ll: iLehtinen et al.l 
l2007l:ISchnee et al.ll2008l) . 

2. As seen in Sect. l2.2[ the NIR extinction map is converted 
into column density using the standard ratio of selective to 
total extinction Ry = 3.1 corresponding to the diffuse ISM. 
However, Ry is expected to increase at high densities (with 
typical extinction Ay > 3), up to abo ut 4.5 owing to grai n 
growth by accretion and coagulation (IWhittet et al.l 1200 iT) . 
This increase lowers the column densities derived from the 
extinction. Therefore the computed values of T25o/Nu could 
be systematically underestimated in dense regions. 

The increase of T250 /Nu by a factor around 2 between the atomic 
phase and the molecular phase obviously tends to decrease the 
equilibrium temperature of the dust particles. For constant in- 
tensity of the incident radiation field and at a first order, the 
temperature of the dust particles should follow the relationship 
T a (t25o/Nu)~^^^^'^^\ which gives for = 1.8 a decrease of 
2K for dust at 17 K. However, the radiation field is attenuated 
in dense regions, and radiative transfer effects must be consid- 
ered for a complete analysis of the spatial variations of the dust 
temperature, which is beyond the scope of this paper. 



6. Conclusions 

Combined with IRAS maps at 100 yum (3000 GHz), HFI maps 
at 857, 545, 353, and 143 GHz allow the precise measurement 
of the emission spectrum of thermal dust with unprecedented 
sensitivity from the faintest atomic regions to the densest parts 
of the Taurus molecular complex. 

While the dust particles located along the lines of sight have 
no reason to be at the same temperature and may have differ- 
ent optical properties, we find that for each pixel of the map the 
measured spectra are reasonably fitted with a single modified 
blackbody, which gives one dust temperature, one spectral emis- 
sivity index, and one dust optical depth per pixel. However, the 
modified blackbody can be slightly broadened around the peak 
of the spectrum because of the range of dust temperature along 
the line of sight, which explains the negative residuals found at 
353 GHz (around -7%). On the other hand, the positive resid- 
uals at 143 GHz (around -hi 3%) could be attributed to a slight 
decrease of the spectral emissivity index at low frequency. 

The dust temperature map we derive from the pixel-by-pixel 
fits provides a spectacular description of the cooling of the ther- 
mal dust across the whole complex from about 17.5 K to about 
13 K. These variations can be caused by variations of both the 
excitation conditions and the optical properties of the dust parti- 
cles. 

The spectral emissivity index map we derive presents signif- 
icant spatial variations, from 1.6 to 2. The distribution is cen- 
tred on 1.78, with systematic error estimated to be 0.07. We 



^ For the same reason, we can also conclude that the inverse correla- 
tion of T250 /Nil with T is robust against the unavoidable combination of 
dust temperature on the line of sight, and reveals a real anti-correlation 
of T250/N11 and T for the dust particles. 
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have checked that the synthetic spectra compute d with the post- 
Spitze r dust models of Draine & Lil (1200 7) and Compiegne et al. 
(|2Q1 1|) have almost identical values of their spectral emissivity 
indexes. Slightly higher values (> 1.8) are found in the coldest 
(about 14 K) structures, and we detect sl T - j3 anti-correlation 
that cannot be explained by the statistical noise and the CIBA. 

We also derive a dust optical depth map with a very high dy- 
namic range, which reveals the spatial distribution of the column 
density of the molecular complex from the densest molecular re- 
gions to the faint diffuse regions. Using the NIR extinction as an 
independant tracer of the column density, we report an increase 
of the measured dust optical depth at 250yum per unit column 
density in the molecular phase by a factor of about 2 compared 
to the value found in the diffuse atomic ISM. The increase of op- 
tical depth per unit column density for the dust particles could be 
even higher in dense regions, owing to radiative transfer effects 
and the increase of Ry in dense regions. 
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Appendix A: 
fitting 



Error propagation - simulation of the 



AA. Principle of the simulations 

We performed Monte-Carlo simulations to understand the 
propagation of the calibration errors, the statistical noise, and 
the CIBA in the determination of the dust temperature T, the 
spectral emissivity index /3, and the dust optical depth T250 at 
250 yum from the fit to a single modified blackbody of the five 
bands at 3000 GHz (100 yum), 857, 545, 353, and 143 GHz 
(Sect. l3.3l) . In practice, we computed 1000 synthetic spectra 
with the same values of T and JS (T = IIK and = 1.8) at the 
central frequency of the five bands. We take a fixed value of the 
dust optical depth (T250 = 1) to study the calibration errors, and 
fixed values of the brightness at lOOyum to study statistical 
noise and CIBA. Then we randomnly add some error for each 
band computed using a Gaussian distribution with a standard 
deviation cr (and mean of 0) equal to the assumed noise. 
Finally, we apply our fitting procedure, and obtain for each 
simulation a set of 1000 values of T,I3, and T250. For the 
diff'erent simulations we conducted, the standard deviations 
obtained for each parameters are given in Table A. 1 , while 
Fig. lA. II shows the correlations between the three parameters. 



A.2. Calibration errors 

In order to study the propagation of calibration errors, we took 
standard deviations of the simulated noise equal to the 
calibration errors presented in Sect. l2.1l For the two HFI bands 
at 545 and 857 GHz, we assumed that the calibration errors are 
fully correlated between the two bands, so in the simulation we 
took the same realization of the simulated Gaussian noise. In 
contrast, the calibration error at 143 GHz is not correlated with 
the calibration error of the two high frequency bands. 
The fitted values of T and JS computed for the 1000 synthetic 
spectra are presented in Fig. lA. II (black symbols). We observe 
the classical anti-corr elation intrin sic to the noise identified by 
several authors (e.g., IShetty et aDl2009 b). The dust optical 
depth at 250/im T250 is also anti-correlated with T owing to the 
temperature dependance of the Planck function. Finally T250 is 
correlated with fi, as expected, since both T-JS and T-T250 are 
anti-correlated. 

To first order, we can consider that for each band the calibration 
erro r is co nstant over the maps, so the errors on T^JS, and T250 
(Fig. lA. II and Table A. 1) systematically aff'ect (in the same 
direction for all pixels) the three parameters derived from the 
fits. Obviously, these calculations are preliminary, since other 
systematic eff'ects are neglected, but they give for this early 
analysis an estimate for the systematic errors on the three 
parameters T,/3, and T250, which are derived from the data. 

A.3. Statistical errors and CIS anisotropies 

The statistical noise on the data is considered to be 
un-correlated both spatially and spectrally. We take standard 
deviations of the simulated noise equal to the statistical noise 
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Fig. A.l. Correlation between the fitted values of T, yS, and T250 
for 1000 simulated spectra at 100 yum 857, 545, and 143 GHz 
(350, 545, and 2100 yum) with T = llK,/3= 1.8, and T250 = 1: 
black, with systematic errors on the gains; red and green, with 
statistical noise and CIBA with 7(100 yum) = 1 and lOMJy sr"\ 
respectively. Relative values are given for T250. 



used for the pixel-per-pixel fitting of the data at 100 yum, 857, 
545, and 143 GHz presented in Sect.O We also take into 
account the noise caused by the CIBA, usin g the standard 
deviations measured by Pla nck/HFl ( Planck C ollaborationI 
i2011n) and IRAS at 100 yum (iPenin et al. 2011 ). 
The correlation diagrams of Fig. lAJl and the standard 
deviations given in Table A. 1 show that the eff'ects of the 
statistical noise and CIBA have a smaller amplitude than the 
eff'ects of the systematic errors, but they depend on the absolute 
brightness, as illustrated in Fig. lA.2l It is interesting to note that 
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Table A.l. Standard deviations of the three parameters T,p and T250 adjusted for different simulations, with T = 17 K and p = 1.8. 



Simulation type 


/(lOOyum) (MJy sr'^) 


o-T (K) 






Calibration errors 




'-0.7 


-0.07 


- 18% 


Statistical noise and CIBA 


1 


~ 1 


-0.25 


-20% 


Statistical noise and CIBA 


10 
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- 0.025 


-2% 
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Fig. A.2. Effects of the statistical noise and GIB anisotropics. Mean values and standard deviations of the fitted values of T, p and 
T250 (relative values) are shown as a function of the 100//m brightness, with T = l7K,j3= 1.8 and T250 = 1. 



above / (100 fim) = 1 MJy sr~^ , the fitted values of T, jS, and 
T250 are not biased, since (T) ^ 17 K, (J3) ^ 1.8, and (t25o> = 1 
(relative value). On the contrary, below 
/(lOOyum) = 1 MJy sr"\ we see that (T) decreases, while (J3) 
and (t25o) increase. Therefore, we consider that pixels with 
I (100 /dm) < 1 MJy sr~^ cannot be used for any quantitive 
analysis of the fitted parameters. 
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